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ABSTRACT 

Dwarf nova outbursts often show coherent (Q ~ 10 4 - 10 6 ) sinusoidal oscillations with the largest pulsed 
fraction in the extreme ultraviolet. Called dwarf nova oscillations (DNOs), they have periods of P ss 3-40 s 
and scale with luminosity as P oc L _/3 with f3 ~ 0. 1 - 0.2. We propose that DNOs may be produced by nonradial 
oscillations in a thin hydrostatic layer of freshly accreted material, the "spreading layer" (SL), at the white dwarf 
(WD) equator. This would naturally explain a number of key properties of DNOs, including their frequency 
range, sinusoidal nature, sensitivity to accretion rate, and why they are only seen during outburst. In support of 
this hypothesis we construct a simple model that treats the SL as a cavity containing shallow surface waves, each 
with the same radial structure, but split into three different modes denoted by their azimuthal wavenumber, m. 
The m = latitudinally propagating mode best matches the periods and scalings associated with most DNOs, 
and DNOs with periods shorter than the WD Keplerian period are explained by the m = -1 prograde mode. 
We also predict a third set of oscillations, produced by the m = 1 retrograde mode, and show its expected 
dependence on accretion rate. 

Subject headings: accretion, accretion disks — novae, cataclysmic variables — stars: oscillations — white 
dwarfs 



1. INTRODUCTION 

Rapid oscillations are often seen in the optical, soft X-ray, 
and extreme ultraviolet (EUV) flux of cataclysmic variables 
(CVs) during dwarf nova (DN) outbursts. These oscillations 
have periods in the range of P as 3 - 40 s, show a monotonic 
relationship with EUV luminosity (Mauche 1996) and thus 
accretion rate, M, and are fairly coherent (Q = \dP/dt | _1 ~ 
10 4 - 10 6 ). They are called dwarf nova oscillations (DNOs; 
Patterson 1981) in contrast to the much less coherent (Q ~ 
5-20) quasi-periodic oscillations (QPOs) also seen from CVs 
and the highly coherent oscillations (Q ~ 10 12 ) that reflect the 
spin of a magnetic primary (e.g., DQ Her). Since their discov- 
ery (Warner & Robinson 1972), DNOs have been seen in ~ 50 
systems (Warner 2004). The DNO periods roughly track the 
surface Keplerian period of the accreting white dwarf (WD) 
[Pk = 2ir/(GM/R 3 y/ 2 ; see Fig. 19 of Patterson 1981, or 
Knigge et al. 1998], so that it is generally believed that the 
oscillations are created near or on the WD surface. 

A number of explanations have been proposed for DNOs, 
but none have quantitatively explained the majority of their 
key features. The characteristic period of the DNOs, P ~ Pk, 
implies a pulsational, rotational, or even perhaps inner ac- 
cretion disk origin. Since the period changes on the same 
timescale as accretion, it is improbable that the entire WD 
takes part in creating DNOs, ruling out global g-modes or 
the WD spin modulated by hot spots. This led Paczyhski 
(1978) to suggest that if accretion torques were causing the 
period drifts then only a small amount of material is involved 
in making DNOs (< 10 24 g). Papaloizou & Pringle (1978) 
studied rotationally modified global g-modes and r-modes on 
WDs. Since truly global modes are not sensitive to the in- 
stantaneous accretion rate, they proposed that DNOs are high 
radial-order modes that have large amplitudes near the WD 
surface. Their work explains many of the properties of DNOs, 
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but some questions still remain, including: (1) what physical 
mechanism favors modes that are strongly concentrated in the 
outer layers of the star, (2) why the radiating area associated 
with DNOs is such a small fraction of the total WD surface 
area (~ 10" 4 to 10" 2 as inferred from the EUV emission of 
SS Cyg during DN outbursts; Cordova et al. 1980; Mauche 
& Robinson 2001, hereafter MR01; Mauche 2004), and (3) 
why the modes are only excited during outbursts. These fail- 
ings led others to suggest alternatives, including magnetic 
accretion onto a slipping belt (Warner & Woudt 2002) and 
non-axisymmetric bulges on the inner part of accretion disks 
(Popham 1999). 

Inogamov & Sunyaev (1999) presented a new way of study- 
ing disk accretion close to an accreting neutron star (NS), 
which Piro & Bildsten (2004; hereafter PB04) extended to 
the case of WDs. Instead of using a boundary layer model, 
they follow the latitudinal flow of accreted material over the 
stellar surface, starting at the equator and streaming up to- 
ward the pole, using a spreading layer (SL) model. In the 
case of WDs the covering fraction (which is equivalent to the 
spreading angle) is found to be small, / ~ 10" 3 to 10 _1 (for 
M w 10 16 - 10 18 g s" 1 ), so that the size of the freshly accreted 
layer is set by the thickness of the accretion disk. 

The material in the SL is much hotter in temperature and 
lower in density than the underlying WD. This contrast allows 
waves in the SL to travel freely, unencumbered by the material 
below. We propose that DNOs are shallow surface waves in 
the layer of recently accreted material confined to the WD 
equator. We argue that the m = mode provides the period 
that is most often identified as a DNO, and we show that the 
m = -1 prograde mode explains the occurrence of some of 
the higher frequency DNOs, including the extra periods seen 
from SS Cyg and VW Hyi (MR01; Woudt & Warner 2002). 
We also discuss the m = 1 retrograde mode and whether it also 
corresponds to an observed frequency. 

In §2 we describe a simple model to estimate the shallow 
surface wave periods and show how they scale with M and 
the WD mass and radius. We then compare these periods to 
observed DNOs from CVs in §3. We conclude in §4 with a 
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summary of the DNO properties successfully understood us- 
ing our model along with a discussion of problems that are 
still left unresolved. 

2. NONRADIAL OSCILLATION PERIODS IN THE 
SPREADING LAYER 

Patterson (1981) showed that (at the time) all known DNOs 
on WDs with measured masses have P > P K - This relation is 
an important constraint for any explanation of DNOs, so we 
begin by showing why a SL mode should mimic such a period. 
The geometry of the SL is a thin, quickly spinning, layer in 
hydrostatic balance, which is predicted both observationally 
(Mauche 2004) and theoretically (PB04) to cover a fraction of 
the surface area / ~ 10" 3 - 10" 1 . In situations where there is a 
large entropy contrast between a surface layer and the under- 
lying material, nonradial oscillations can be confined to high 
altitude regions with little or no pulsational energy extending 
deeper into the WD. When the horizontal wavelength is much 
greater than the layer depth, these shallow surface waves have 
a frequency 
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(1) 

where g eff = GM/R 2 -v\jR- v 2 /R w GM/R 2 - v\jR is the 
surface gravitational acceleration decreased by centrifugal ef- 
fects, h = P I (pgeff) is the pressure scale height of the layer, and 
k is the transverse wavenumber. The SL can be thought of as 
a waveguide with latitudinal width 2fR and azimuthal length 
2-jtR, but since / <C 1 the latitudinal contribution dominates 
so that k ~ l/(fR). Setting g eB = XGM/R 2 , where A < 1 is a 
dimensionless parameter that depends on the spin of the layer, 
we rearrange the terms in equation (1) to find 
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which shows that the mode's frequency is the Keplerian fre- 
quency times a factor less than unity (as long as Xh < f 2 R), 
and therefore explains why these shallow waves are consistent 
with the findings of Patterson (1981). 

To find how these modes scale with M and the WD mass we 
consider a simple model to estimate / and h. The theoretical 
studies of PB04 suggest that in the M range of DN the cov- 
ering fraction, /, may be set by the thickness of the accretion 
disk resulting in (Shakura & Sunyaev 1973; in the limit of 
gas pressure much greater than radiation pressure and using 
Kramer's opacity) 



/=1.8xl0- 2 a^ 2 X /2 X 3/8 < 8 ; 



(3) 



where oyisk is the viscosity parameter for the accretion disk 
and a dlsk , 2 = M xl =M/(10 17 g s" 1 ), M\ =M/M Q , 

R 9 = R/(10 9 cm), and we set the factor [1 -(r/PT 1 / 2 ] w 1. 

As long as the WD is rotating at much less than breakup 
(Gansicke et al. 2001), half of the accretion luminosity is 
released in a layer at the WD surface, so that the flux of 
this layer is given by AirfR 2 F = GMM/(2R). In a one-zone 
layer the radiative flux equation can be integrated to give 
F = acT 4 /(3ny), where a is the radiation constant, k is the 
opacity (set to 0.34 cm 2 g" 1 for Thomson scattering with a 
solar composition), and y is the column depth of the layer 
(measured in units of g cm" 2 ). The column depth is set by 
continuity to be y = Mf v i sc /(47r/P 2 ), where f v ; sc = h 2 jv is the 
timescale for viscous dissipation in the SL. The viscosity be- 
tween the fresh, quickly spinning material and the underlying 



WD is given by v = asiy<f>h (PB04), where asL is the SL's vis- 
cosity parameter, and w (GM/R) 1 / 2 is the initial azimuthal 
velocity of the material in the SL. Using an ideal gas equa- 
tion of state, we find the lowest order mode, which does not 
propagate in the azimuthal direction so we denote it with an 
azimuthal wavenumber m = 0, has a period 
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where asL,3 = o:sl/10~ 3 and Ai = A/10" 1 . Equation (4) pro- 
vides both a scaling and period suggestive of DNOs. We dis- 
cuss these similarities in more detail when we compare to ob- 
servations in §3. If we instead take v = ots L c s h (in analogy to 
a Shakura & Sunyaev accretion disk), where c s is the speed 
of sound, we find P oc m- 13 / 140 M- 27 / 56 P 89 / 56 , also similar t0 
observations. This shows the robustness of this idea in repli- 
cating the general scalings of DNOs. 

The arguments from above suggest that the SL can con- 
tain additional modes, most notably those that propagate in 
the azimuthal direction. These modes have an observed fre- 
quency of w bs = w-oto;sl|, where wsl is the spin of the SL 
and m is the azimuthal wavenumber. Since the layer is spin- 
ning quickly, it is possible that Coriolis effects will modify 
k, which has been studied in the limit of a thin layer using 
the "traditional approximation" (see Bildsten, Ushomirsky & 
Cutler 1996 and references therein). Using this analysis, Cori- 
olis effects are negligible since / < cj/(2wsl)- If we set 
Psl = 27t/wsl, where Psl ^ Pk, the next lowest order modes 
(still n = 1) have periods 
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Since there are a few DNOs with P < Pk (as we discuss in 
§3), we propose that the prograde m = -1 mode is also neces- 
sary for explaining some short period DNOs. It is less clear 
whether the m = 1 retrograde mode is consistent with an oscil- 
lation observed during a DN outburst, and in §3 we speculate 
whether this mode may be related to the long-period DNOs 
(lpDNOs; Warner, Woudt & Pretorius 2003). 

3. DNOS IN THE CV POPULATION 

We compare the most recent compilations of DNO peri- 
ods (Table 1 of Warner 2004) and WD masses (Ritter & Kolb 
2003) of CVs in Figure 1 . The period error bars indicate the 
range of periods that have been seen from each object, while 
the mass error bars indicate the errors. We plot systems that 
have shown both high and low DNO periods as two separate 
points (these CVs are SS Cyg, CN Ori, and VW Hyi). The 
thick, dashed line denotes the surface Keplerian period, Pk, 
for a given WD mass (using the mass-radius relation of Truran 
& Livio 1986). This demonstrates that there are DNOs both 
above and below P K , which is why we consider SL nonradial 
oscillations with m = and m = -1 . For each mode we present 
periods for accretion rates M = 10 16 - 10 18 g s" 1 , the range ex- 
pected during a DN outburst. To calculate the m = —\ mode 
we must assume a period for the SL's spin, Psl, so we plot 
mode periods for both P SL = Pk and P SL = 2P K . This shows 
that our model is insensitive to the choice of the SL spin rate. 

Besides predicting multiple classes of DNOs with different 
period ranges, we also predict that each of these groups will 
have a different dependence on M. This can be seen from 
the size of each shaded band, which is much wider for the 
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FIG. 1. — DNO periods vs. measured WD masses. The vertical bars show 
the range of periods observed from each CV, while the horizontal bars cor- 
respond to the mass measurement errors. CVs that show two distinct DNO 
periods are plotted as two separate points. The thick dashed line is the Kep- 
lerian period, Pk. The m = mode (wide, dark-shaded region) is plotted for 
M = 10 16 — 10 18 g s" 1 (from top to bottom), using eq. (4). The m = —1 mode is 
also plotted for M = 10 16 - 1() 18 g s _1 (from top to bottom), for both P SL = P K 
(light-shaded region) and Psl = 2Pk (medium-shaded region), using eq. (5). 
The m = -1 modes show less variation with M in comparison to the m = 
mode. 



m = mode than the m = -1 mode. On average the DNOs with 
P < Pk show less variation with M, qualitatively consistent 
with the m = — 1 mode of our model. 

To investigate the M dependence in more detail we plot our 
predicted DNO periods versus M for an M = 1 .OM WD in 
Figure 2. This illustrates the shallower dependence for the 
m = — 1 mode. The separation and relative slopes of the m = —l 
and m = modes are suggestive of the frequency doubling that 
was seen by MR01 during the 1996 October outburst of SS 
Cyg. In both the beginning and tail of the burst, they saw a 
DNO with P = 6.59 - 8.23 s and following a P oc IT? power 
law with [3 = 0.097. Near the burst peak, P suddenly shifted 
to 2.91 s and then decreased with a shallower power law with 
(3 = 0.021. Indeed, we find similar periods and power laws 
when we use a larger WD mass. 

In Figure 2, we also plot the m = 1 retrograde mode for 
Psl = (1 — 1.5) x P K . Owing to the difference taken in equa- 
tion (5), this mode has a much more complicated dependence 
on M, no longer being a power law nor monotonic. This mode 
has a higher period than typical DNOs, and may be relevant 
for the lpDNOs. Warner (2004) identified 17 CVs as showing 
such oscillations. The lpDNOs typically have larger ampli- 
tudes than DNOs and, like DNOs, are not seen in every DN 
outburst. Often, both of these oscillations are seen simultane- 
ously. To positively identify these as modes requires checking 
for the predicted scalings between P and M. This test may 
negate a mode explanation for lpDNOs since Warner, Woudt 
& Pretorius (2003) claim to find no such correlation for these 
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FIG. 2. — SL nonradial oscillation period compared to M for an M = l.OM 
WD. The thick solid line is the m = mode, while the shaded region shows 
the m = -1 mode for a range of Psl = (1.0-1.5) X Pk (bottom to top). The 
dashed lines show the m = 1 mode for the same range of P SL . 



oscillations and thus favor a spin-related mechanism. On the 
other hand, on separate occasions there have been 32 - 36 s 
(Robinson & Nather 1979) and 83- 110 s (Mauche 2002b) 
oscillations seen from SS Cyg, in the domain expected for 
these modes. This wide spread of periods may be explained 
by the m = 1 mode's steeper dependence on M. The 32- 36 s 
oscillations showed much less coherence than typical DNOs, 
which could be a result of the retrograde oscillations beating 
against the accretion disk. Other CVs such as VW Hyi show 
similar modes in the range of ~ 90 s (Haefner, Schoembs & 
Vogt 1977, 1979), which may be of similar origin. 

4. DISCUSSION AND CONCLUSIONS 

We propose that DNOs in outbursting CVs are nonradial 
oscillations in a hot layer of freshly accreted material near the 
WD equator. A large number of DNO properties are then sim- 
ply understood: (1) the highly sinusoidal nature of the oscilla- 
tions is consistent with nonradial oscillations, (2) the periods 
can change on the timescale of accretion because there is lit- 
tle mass in the layer (< 10 21 g; PB04), (3) the periods vary 
inversely with M because they have the temperature scaling 
of shallow surface waves, (4) the covering fraction is natu- 
rally small for the SL, (5) the DNOs are only seen during 
DN outbursts because this is the only time when an optically 
thick layer of material can build up at the equator, and (6) the 
largest pulsed amplitude is in the EUV, consistent with the SL 
temperature. 

In support of our hypothesis, we presented a simple, phe- 
nomenologically motivated model that quantitatively explains 
many of the DNO's features. We compared our model with 
the overall population of DNOs, both as a function of WD 
mass and accretion rate. In each case, it is necessary to con- 
sider both the m = mode and the m = -1 (prograde) mode 
to explain the range of periods observed and the different P- 
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M scalings. The majority of DNOs are consistent with the 
period of the m = mode, and this may be related to its lat- 
itudinal propagation. The speed of shallow surface waves is 
ps gh so that they slow down as the scale height decreases. 
This steers an initially azimuthally propagating mode to in- 
stead travel perpendicular to the SL edge, similar to ocean 
waves at the beach, and may explain why m = modes are fa- 
vored over m ^0 modes. In the 1996 DN outburst of SS Cyg 
(MR01), the m = -1 mode is only seen near the outburst peak. 
This may indicate that a wide enough region without differen- 
tial rotation is only present when there is sufficient spreading 
at high M. 

There are still many difficulties that must be addressed 
about this idea for explaining DNOs. From the SL model 
(PB04) we borrowed the concept of hot material in hydro- 
static balance covering a small fractional area of the WD, but 
we ignored important details of these calculations, such as 
differential rotation and the change of scale height with lati- 
tude. The data do not require such additional complications, 
so we refrain from including them for now. In a more sophis- 
ticated model, it would probably still be true that kocl/ ( fR), 
but an eigenvalue calculation would determine the constant of 
proportionality. A full calculation may also help explain the 
variety of power laws, /3 = -dlogP / dlogL, for various CVs, 
since our model only predicts one unique power law index for 
the m = mode. The WD spin (which has so far been ne- 
glected) may be an important differentiating variable as well. 

Another important property of WDs that would affect the 
modes is a magnetic field. A strong field inhibits shearing be- 
tween the SL and WD and modifies the frequency of shallow 
surface waves. It is therefore interesting that no intermedi- 
ate polars (IPs) have shown DNOs or lpDNOs. Even LS Peg 
and V795 Her, neither of which are IPs, but both showing 
polarization modulations (Rodriquez-Gil et al. 2001, 2002) 
indicative of a reasonably strong magnetic field, are with- 
out DNOs or lpDNOs. Further studies should also work to- 



ward an understanding of the excitation mechanism for the 
modes, which would help explain the high coherence typical 
of DNOs. Material deposited at the WD equator spreads fairly 
quickly, 10 - 100 s (PB04), so that this material must be mov- 
ing through the oscillating region on timescales of order the 
mode period. In light of the many strengths of the SL mode 
explanation of DNOs we do not abandon it because of this 
problem, but this is definitely a weakness that puts limits any 
proposed excitation mechanism. 

Our explanation of DNOs raises interesting questions about 
the relationship between oscillations originating from accret- 
ing compact objects. Mauche (2002b) showed that there is a 
strong correlation in the high to low oscillation frequency ra- 
tio of WDs, NSs, and black holes (BHs). Using this picture, 
DNOs are associated with the kilohertz QPOs of low mass 
X-ray binaries. Suggestively, the Fourier frequency resolved 
spectroscopy of NSs (Gilfanov, Revnivtsev & Molkov 2003) 
imply that both the normal branch oscillations and the kilo- 
hertz QPOs are created in the NS boundary layer, strength- 
ening the possibility that our SL mode model may also apply 
in this case. For a typical NS mass and radius our model re- 
sults in frequencies in the range expected for kHz QPOs, but 
it does not explain their complicated P-L relation (the "paral- 
lel tracks;" van der Klis 2000). Interestingly, DNOs may also 
show the parallel tracks phenomenon, as seen in three obser- 
vations of SS Cyg (Mauche 2002a), once again supporting the 
correlation. On the other hand, continuing the analogy to WD 
and BH oscillations is problematic because in the case of BHs 
there is no surface for nonradial oscillations. 
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